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Abstract. We study interacting systems of linear Brownian motions whose 
drift vector at every time point is determined by the relative ranks of the coor- 
dinate processes at that time. Our main objective has been to study the long 
range behavior of the spacings between the Brownian motions arranged in in- 
creasing order. For finitely many Brownian motions interacting in this manner, 
we characterize drifts for which the family of laws of the vector of spacings is 
tight, and show its convergence to a unique stationary joint distribution given 
by independent exponential distributions with varying means. We also study 
one particular countably infinite system, where only the minimum Brownian 
particle gets a constant upward drift, and prove that independent and identi- 
cally distributed exponential spacings remain stationary under the dynamics 
of such a process. Some related conjectures in this direction have also been 
discussed. 

1. Introduction 

In this paper we consider systems of interacting one-dimensional Brownian mo- 
tions X = (Xi(t), i E I, t > 0), where i ranges over an index /, which is either the 
finite set {1, . . . , N}, or the countable set of positive integers N. 

For I = {1, 2, ... , N}, if we define ordered coordinates of any vector x £ M. N by 

X(l) < X( 2 ) < • • ■ < X(JV), 

the locations Xi(t) of the Brownian particles evolve according to the system of 
stochastic differential equations 

(1) dXi(t) = ^Sjl { Xi{t) = x u) {t) } dt + dB,{t) (i e I) 

for some sequence of drifts Si, §2, ■ ■ ■ S R. Here the Bi's are asumed to be inde- 
pendent (.F^-Brownian motions for some suitable underlying filtration (Tt). Less 
formally, the Brownian particles evolve independently except that the zth ranked 
particle is given drift Si. 

For these finite systems, with arbitrary initial distribution of (Xj(0), i € J), and 
arbitrary drifts Si, the existence and uniqueness in law of such an N particle system 
is guaranteed by the standard theory of SDE's. See Lemma[6]for more details. We 
would also like to define SDE (fT]) when I is the entire countably infinite collection 
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N. This is trickier since the ordered process might not remain well-defined anymore, 
and the existence of the solution of SDE ([I]) depends both on the initial distribution 
and the sequence of drifts. We consider one such system in Section [3] with a drift 
sequence (<5, 0, 0, . . .) (5 > 0) for which the finite-dimensional arguments can be 
suitably extended. For this system to exist in the weak sense, starting from an 
initial vector Xi(Q) < ^(0) < . . ., we show that it suffices to assume 

liminf Wo)-x l( o)) 2>o 

See Lemma [TT] for the details of the proof. 

Consider now the Brownian spacings system derived from this ordered Brownian 
particle system with rank dependent drifts, i.e., 

A k (t) := X {k+l) (t) - X {k) (t) for k, k + 1 G /. 

The ordered particle system derived from independent Brownian motions with 
no drift (meaning Si = 0) was studied by Harris [15j . Arratia [T] and Sznitman 
[IS] j [301 P- 187]. By Donsker's theorem, this system can be interpreted as a scal- 
ing limit of ordered particle systems derived from independent symmetric nearest 
neighbour random walks on Z. Harris |15j considers the spacings of an infinite 
ordered Brownian particle system defined by : 

A*(t) := B {t+1) (t) - B {t) (t) (ieZ) 

where {Bi} is a family of independent Brownian motions with no drifts and initial 
states .Bj(O) = S(;)(0) such that B (0) = and the B 4 (0) for i G Z\{0} are points 
of a Poisson process of rate A on R. That is, (-B(j) (t),t > 0) ^ is the almost surely 
unique collection of processes with continous paths such that B^ (t) < -B(i+i) (t) 
for all i G Z, t > and the union of the graphs of these processes is identical 
to the union of the graphs of the Brownian paths (Bi(t),t > 0). ^. We call 

(Bu\(t),i eZ,(>0) the Harris system of ordered Brownian motions with rate A 
and their differences (A*(t),i eZ,t>0) the Harris system of Brownian spacings 
with rate A. 

As observed by Arratia [5J §4], the corresponding stationary system of spac- 
ings between particles of the exclusion process associated with a nearest-neighbor 
random walk on Z can be interpreted as a finite or infinite series of queues, also 
known as the zero range process with constant rate. See [THJ [17l \TE\ [20l [T9l [24] for 
background on systems of Brownian queues. Such connections between systems of 
queues and one-dimensional interacting particle systems have been exploited by a 
number of authors, in particular Kipnis [35], Srinivasan [25], Ferrari- Fontes [13] . 
[14] . and Seppalainen [57]. Ferrari [T^] surveys old and new results on the limiting 
behavior of a tagged particle in various interacting particle systems. Also see arti- 
cles by Baryshnikov [4] and O'Connell-Yor [24] for some recent studies of Brownian 
queues in tandem, connected to the directed percolation and the directed polymer 
models, and the GUE random matrix ensemble. 

More recently rank dependent SDEs have been considered by several authors 
as possible models for financial or economic data. Fernholz in ID] introduces the 
so-called Atlas model which we study in this paper. It is a model of finitely many 
Brownian particles, where at every time point the minimum Brownian motion gets 
a constant positive drift, while the rest of them get no drift. The general rank 
dependent interacting Brownian models, whose drifts and volatilities depend on 
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time- varying ranks have been considered by Banner, Fernholz, and Karatzas in [3J, 
with whom our work in this paper bears close resemblance. For SDE ([1]) , they work 
under a specific condition on the drift sequence required for stability of the solution 
process. We prove in Theorem [8] (condition (jT3j) ) that this condition is indeed 
necessary and sufficient. Although their method is mostly based on a beautiful 
analysis of the local times of intersections of different Brownian motions, they 
also note the connection with the Harrison- Williams theory of reflected Brownian 
motions which we use extensively in this paper. See Lemma 2] for a complete 
statement. 

In that same article [3J the authors establish marginal convergence of spacings 
Afe to exponential distributions. They leave the question of joint convergence open, 
which we settle in this paper in Theorem [8] by proving that the vector of spacings 
converge jointly to independent exponentials with different means. They also study 
ergodic properties of such processes including a demonstration of the exchange- 
ability of the indices of the Brownian particles under rank-dependent drifts and 
volatilities. In their later papers, Fernholz and Karatzas also consider generalized 
versions of where the drifts and the volatilities depend on both the index and 
the rank of a Brownian particle. As expected, in most such cases, explicit descrip- 
tions of their properties become very difficult to obtain. However, many interesting 
results can still be recovered. A good source for what has been done so far can 
be found in the recent survey article by Fernholz and Karatzas [11]. Also see the 
article by Chatterjee and Pal [5] for a follow-up in this direction where the authors 
consider an increasing number of Brownian particles in a rank-dependent motion 
and establish connection with the Poisson-Dirichlct family of point processes. 

McKean and Shepp, in [23j . also consider Brownian motions interacting via their 
ranks. They start with two Brownian motions, and their objective is to find the 
optimal drift under constraints (as a control) such that the probability that both 
Brownian motions never hit zero is maximized. The solutions is, as they establish, 
the Atlas model for the two particle system. 

An interesting related model, studied by Rost and Vares in replaces the 
ordered particles in the Harris model by linear Brownian motions repelled by their 
nearest neighbors through a potential. The authors study stationary measures for 
the spacings of such processes, and show that rescaled combinations of spacings 
converge to an Ornstein-Uhlenbeck process. 

Our purpose here is to draw attention to the general class of Brownian particle 
systems with rank dependent drifts, as considered in [3J. Many natural questions 
about these systems remain open. We are particularly interested in an infinite 
version of the Atlas model, with a drift sequence (S, 0, 0, ...) with Atlas drift 6 > 0. 
One result we obtain for this system is the following: 

Theorem 1. For each S > 0, a sequence of independent Exponential(2S) variables 
provides an equilibrium distribution of spacings for the infinite Atlas model with the 
Atlas drift 5. 

Theorem [TJ which essentially follows from Theorem [TU suggests a number of 
interesting conjectures and open problems. In particular, we can immediately for- 
mulate: 

Conjecture 2. For each S > 0, Theorem^ describes the unique equilibrium distri- 
bution of spacings for the infinite Atlas particle system with Atlas drift S. 
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Let (Ai(t), A2(t), . . .)t>o denote the equilibrium state of spacings of the infinite 
Atlas Brownian particle system described by Theorem [TJ This process has some 
subtle features. For each k — 1,2, . . . and each t > 0, 

(2) (A fe (i),A fe+1 (t),...) i (A 1 (t),A 2 (t),...) 

and the common distribution of these sequences is that of independent exponen- 
tial (26) variables. But while both sides of ([2]) define stationary sequence-valued 
processes as t varies, these processes do not have the same law for all k. In 
particular, the finite-dimensional distributions of non-negative stationary process 
(Afc(i), t > 0) depend on k. 

Harris [T51 eqn. (7.1)] gave an explicit formula for the law of (t), the location 
at time t of particle initially at in the Harris system of ordered Brownian motions, 
from which he deduced for 2A = 1 that 

(*\ £(Q)W d (2\ 1/4 1 

(3) -M^u) vrx {) " 



where B(l) is standard Gaussian. As remarked by Arratia [TJ p. 71], the conclusion 
for general A > folows from the case 2A = 1 by Brownian scaling. See also De 
Masi, Ferrari [6], Rost, Vares [26], and Arratia [TJ where variants (or generalizations) 
of is proved for a tagged particle in the exclusion process on Z associated with 
a simple symmetric random walk. Harris conjectured that the process i?( ) is not 
Markov, and left open the problem of describing the long-run behaviour of paths of 
_B( ) ■ These questions were answered by Diirr, Goldstein, and Lebowitz in [71 Thm 
7.1] where they prove that the tagged process (t), suitably rescaled, converges to 
fractional Brownian motion with Hurst parameter 1/4. In fact, they show a general 
theorem where convergence to fractional Brownian motion holds for systems with 
processes which have stationary increments and perform elastic collisions as in the 
Harris model. Also see [8] where the same authors generalize their results in the 
case where an external potential is present. 

In this connection, we have the following conjecture: 

Conjecture 3. For each 6 > 0, in the infinite Atlas model with Atlas drift S > 0, 
and initially ordered values X^(0), k — 1, 2, . . . which are the points of a Poisson 
process with rate 26 on (0,oo), 

A (fc) (t)-A (fc) (0) d f2\ 1/i c k 

tTa * — B(l) as t — > oo 

i 1/4 W V26 

for some sequence of constants c k > with — > 1 as k — > oo. 

The paper is organized as follows. In the next section, we analyze the finite 
Brownian particle system with rank dependent drifts. The main result is Theorem 
[HI which describes the convergence in total variation of the laws of the spacings 
to that of independent exponential distributions. The precise condition needed on 
the drift sequence for such stability has also been proved. In Section [3] we look at 
countably infinite Brownian particles with the dynamics of the Atlas model. The 
main result, Theorem [TJ follows readily from Theorem 1141 

2. The finite Brownian particle system 

We first present some results regarding the asymptotic behavior of spacings for 
the N particle system with arbitrary rank-dependent drifts 6i,l < i < N. We 



BROWNIAN MOTIONS INTERACTING THROUGH RANKS 



5 



start by recording the following two lemmas, which clarify the issues of existence 
and uniqueness of the N particle system with arbitrary drifts Si, and character- 
ize the associated ordered particle system. See the book by Revuz and Yor [55] 
for a background and the definitions of concepts from the calculus of continuous 
semimartingales . 

Lemma 4. Let Xi,l < i < N be a solution of the SDE (Qp defined on some filtered 
probability space (Q,J 7 ,{J-'t},P) ) for some arbitrary initial condition and arbitrary 
drifts {Sj}. Then for each 1 < j < N the jth ordered process Xu\ is a continuous 
semimartingale with decomposition 

(4) dX {j) (t) = dfy(t) + -±=(dL {j _ ltj) (t) - dL UJ+1) (t)) 

where the f3j 's for 1 < j ' < N are independent {Tt) -Brownian motions with unit 
variance coefficient and drift coefficients Sj, and where £(o,i) = -^(Ar.jv+i) — 0, and 
for 1 < j < N - 1 

(5) L 0w) (*)=liml^ l((X {j+1) (s)-X U) (s))/V2<e)ds (t > 0) 

which is half the continuous increasing local time process at of the semimartingale 
(Xu+\\ — Xtj\)/y/2. Moreover, the ordered system is a Brownian motion in the 
domain 

(6) < i<N)£R N : x (1) < x (2) <■■■< x {N) } 

with constant drift vector (5j ,1 < j < N) and normal reflection at each of the N — 1 
boundary hyperplanes {xu\ = for 1 < i < N — 1. 

Proof. Sznitman [551 page 594], [Ml Lemma 3.7] gave these results in the case of 
zero drifts. As he observed, they follow from Tanaka's formula [25; page 223] and 
the definition of Brownian motion in a polyhedron with normal reflection. See, for 
example, Varadhan and Williams |31j . Sznitman further proves that the Brownian 
motions {f3j} are independent, which is essentially due to Knight's theorem |25[ 
page 183]. The factors of \[2 are most easily easily checked in the case N = 2. 
That they are the same for all N follows by a localization argument. The results 
of the lemma for general drifts Sj are deduced from the results with no drift by 
application of the next two lemmas. □ 

Lemma 5. Let X\, X%, ... be a sequence of independent Brownian motions indexed 
by a finite set I . Consider the ordered process Xn\, Xp), ■ ■ ■■ Then, the following 
holds almost surely: 

M ,t 

(7) J2 / l{X i (s)=X {j) (s) }ds = t, Vj < + Vte [0,oo). 

Moreover, the points of increments of the finite variation processes L(j j+i) are 
almost surely disjoint. 

When I is countable, equation (JT]) still holds true as long as the assumptions on 
the initial values ^i(O) < Xa(0) < . . . are sufficient to guarantee the existence of 
the ordered processes for all finite times. 
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Proof. For any two indices k < I, the lebesgue measure of the set { t € [0, oo) : 
Xk = X\ } is zero. This follows because the zero set of Brownian motion is of 
lebesgue measure zero almost surely. 

Now, the event Ylili Jo •"■ { -^«( s ) = ^-(j)( s ) } ds > t for some j and some t 
implies that there exists some pair (k, I) such that the lebesgue measure of the set 
{0 < s < t : Xk(s) ~ X(j)(s) = Xi(s) } is positive. This is of measure zero 
according to the previous paragraph and by countable additivity. 

The other possibility of Y^ili Jo { ^»( s ) = ^(j)( s ) } ds < t is trivially ruled 
out by our assumption that the ordered processes are always achieved. 

For the second assertion, note that, according to the general theory of semi- 
martingale local times [25] , the process Lq j+1 ) increases only on the random closed 
set of times t when X^(t) = Xu + x\(t). These random sets are almost surely 
disjoint as j varies. This is because, with probability 1, there is no triple collision, 
meaning a time t > at which Xi(t) = Xj(t) = Xk(t). It follows from the fact that 
the bivariate process {Xi — Xj, Xj — Xk) is a linear transformation of a standard 
planar Brownian motion which does not hit points. □ 



Call (Q, T, { •?"*}, P) the canonical setup if (Q, T , {•?"*}, P) is the usual space of 
continuous paths in R^, with the usual right-continuous filtration, and Xi is the 
ith coordinate process. 

Lemma 6. Let 5 = (Si, 1 < i < N) € R N and let [i be a probability distribution on 
R N . 

(i) In the canonical setup there is a unique probability measure P' 5 '' 1 under which the 
coordinate processes (Xi, 1 < i < N) solve the system of SDE's ([1]) with initial dis- 
tribution pi. In particular, for <5 = 0, the law P 0,AI is the Wiener measure governing 
standard Brownian motion in M. N with initial distribution [i. 

(ii) In the canonical setup, for each t > 0, the law P S '^ is absolutely continuous with 
respect to P '' 1 on Tt, with density 

(N N 

where which is the same as in can be defined by the expression 

(9) / l{X i ( S )=X (j) (s)}dXi(s), l<j<N. 

i=i Jo 

Under P 5,A1 the (3j 's are independent Brownian motions on K with drift coeffients 
{Sj} and unit diffusion coefficients. 

(iii) // (Xi, 1 < i < A^) is a realization of the N particle system with drifts {Si} and 
initial distribution /i on an arbitrary probability space (fl,J-,¥), then the P joint 
distribution of the processes Xi is identical to the P S '^ distribution of the coordinate 
processes on the canonical space, as specified in (ii). 

Proof. This is an instance of a well known general construction of the solution of 
an SDE with drift terms from one with no drift terms [25l Chapter IX, Theorem 
(1.10)]. 

Under P ^, the fact that {f3j} is a collection of independent Brownian motions 
also follows from expression ((9]) and equation ((7|). Under P the process (3j is a 
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stochastic integral with respect to Brownian motions with drifts. By expanding the 
drift term, it is obvious that (3j is a Brownian motion itself with drift 5j . □ 



Note that the SDE defining the N particle system with drifts is a typical example 
of an SDE for which there is uniqueness in law, but not pathwise uniqueness. 
For a solution X of SDE U), let 

N N 

*(*> = *£*(*>(*) 

1=1 j=i 

which is the center of mass of the particle system, where we regard each particle 
as having mass 1/N. We call the R^-valued process 

(Xi-X,l<i<N) 

the centered system. Note that the N — 1 spacings defined by taking differences of 
order statistics of the original system are identical to the N — 1 spacings defined by 
differences of the order statistics of the centered system. So the N order statistics 
of the centered system, which are constrained to have average 0, are an invertible 
linear transformation of the N — 1 spacings of original system. 

Lemma 7. For the N particle system, with arbitrary drifts and initial distribution, 
(i) the center of mass process is a Brownian motion with drift 

N 

(10) 8 N := N-^Sj 

j'=i 

and diffusion coefficient 1/N; explicitly 

(11) X(t)=X(0)+5 N t+^ 

V N 

1 N 

where B(t) :— N 2 y\ —1 Bj(t) is a standard BM on R; consequently 
X(t) 



t 



Sn almost surely as t — > oo. 



(ii) The shifted center of mass process X(t) — X(0) is independent of the centered 
system ( Xj(-) - X{-) ),l<i<n. 

Proof. Part (i) follows plainly from the SDE |1| by summing over i. 

For part (ii), note that if o~\ and o~i are two independent sub-cr-algebras of a 
probability space (A, A, P), and we change P to another probability Q via defining 
dQ/dP — fg, where / G <7i, g G 02, then o\ and 02 remain independent under Q. 

Now as in Lemma |6l when (Xi) denotes N independent Brownian motions with 
initial distribution /i, the shifted average X(t) — X(0) and the centered process 
Y = X — XI are independent. This follows from the fact that conditionally on 
X , the processes are Gaussian with zero covariance, and that X(t) — X(0) is 
independent of X . Now to get to P 5 ^, the Radon-Nikodym derivative is given by 
dHJ. Note that, from expression © we get 

= E / 1 = y 0)( fi )} dY t (s)+X(t) - X(0). 

i=i ^° 
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Thus, from jHJ, it is clear that dP s ^/dP ^ can be written as fg, where / £ a(Y) 
and g £ a(X(-) — X(Q)). Now, by first localizing at finite time, and using the 
argument in the preceding paragraph, we establish the claim in part (ii). □ 



Theorem 8. For 1 < k < N let 

k 

(12) a k :=J2( S i~^) 



i=l 

where 5n is the average drift as in (|10[) . For each fixed initial distribution of the N 
particle system with drifts {Si, i = 1,2,..., N}, the collection of laws of X( N ^ (t) — 
Xru (t) for t > is tight if and only if 

(13) a k > for all 1 < k < N- 1, 

in which case the following results all hold: 

(i) The distribution of the spacings system (Xy + i) — Xu\, 1 < j < N — 1) at time t 
converges in total variation norm as t — > oo to a unique stationary distribution for 
the spacings system, which is that of independent exponential variables Yj with rates 
2a j, 1 < j < N — 1. Moreover, the spacings system is reversible at equilibrium. 

(ii) The distribution of the centered system at time t converges in total variation 
norm as t — > oo to a unique stationary distribution for the centered system, which 
is the distribution of 

(5,(0-1 - S, 1 < % < N) 
where Sq := and Si := Y\ + • • • + Yi for 1 < % < N — 1, where tt a uniform random 
permutation of {1, ... , N} which is independent of the Yi, and 

N JV-1 JV-1 

i—l i—1 i—1 

Moreover, the centered system is reversible at equilibrium. 

(iii) As t — > oo 

Xi(t)/t — > Sn almost surely for each 1 < i < N, 
and the same is true for X^{t)/t instead of Xi(t)/t. 

Remark. Regard the system as split into left-hand particles of rank 1 to k and right- 
hand particles of rank k+1 to N. If these two parts of the system are started at some 
strictly positive distance from each other, they evolve independently like copies of 
the fc-particle system and the (N ~ fc)-particle system respectively, until the first 
time there is a collision between a left-hand particle and a right-hand particle. It 
follows by summing up over the corresponding drifts that the centers of mass of 
the two subsystems left to themselves would have almost sure asymptotic speeds 
5k and 5k respectively, where 

k n 

4 -.= -r s i and ^ - ( N - fc r x E 

1=1 i=k+l 

Since 

j, r ,c k(N-k) f - % 
ak = kdk - kd N = — [d k - d k ) 
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we see that ctk > iff 5k > 5k, which ensures that the righthand system cannot avoid 
an eventual collision with the lefthand system. According to (fT2|) . for arbitrary 
prescribed Sn G R, and otk > 0, the unique drift vector determining an ergodic 
TV-particle system whose average drift is 5n and whose asymptotic spacings are 
independent exponential variables with rates 2ak, is given by 

(14) Si =6 N + ai ~ cti-x (l<i<N) 

where ao ■— oln '■— 0. Given an arbitrary cumulative probability distribution 
function F on the line, and arbitrary 5 £ R and e > 0, it is clear that by taking 
N suitably large, we can choose (at, 1 < k < N — 1) and hence [5j, 1 < j < N) so 
that for all sufficiently large t 



( 1 N 

V x i=i 



F(x) 



> £ < £ 



and X(t)/t — > 6 almost surely as t — > oo. Thus no matter what its initial distribu- 
tion, the N particle system looks asymptotically like a cloud of particles with mass 
distribution close to F, which is drifting along the line at speed 5. 

We would also like to mention here that a special case of the above theorem has 
been considered in a recent article by Jourdain and Malrieu [21j . They consider 
SDE ([T]) with an increasing sequence of Si's and establishe joint convergence of the 
spacing system to independent exponentials as t goes to infinity. 
Proof. [Proof of Theorem[8] According to LemmalU the ordered N particle system 
is a Brownian motion (X^)(t),l < k < N)t>o in the domain ([6]) with identity 
covariance matrix, constant drift vector (5j , 1 < j < N) and normal reflections 
at each of the N — 1 boundary hyperplanes {x^ = x^ + i^} for 1 < i < N — 1. 
Note that the vector (1, 1, . . . , 1) of all ones is in the intersection of these boundary 
hyperplanes. 

Let (3 be an independent one-dimensional Brownian motion ((3q = 1) with a 
negative drift —9 (6 > 0) which is reflected at the origin. The process defined by 

Y({t) := Xi{t) - X(t) + J=, i = 1,2, ... N, 
v Jy 



has the same spacings as the original process X. Moreover, by the independence of 
the centered system and the center of mass process established in Lemma [71 and 
the fact that Y' — fl/VN is the projection of Y' on the subspace generated by the 
vector (1, 1, . . . , 1), it follows that the vector (Y[, Y^ , . . . , Y' N ") is a Brownian motion 
with identity covariance matrix which is normally reflected in the wedge 

( N \ 



y e R N : J2 Vi ^ °' f i ^ 2/2 < ■ • ■ < VN 



The drift vector in this wedge can be written as a linear combination of spacings 
by summation by parts: 



N / n x N-l / k \ n N 



(15) E [Si - S N - _ j yi = - g ^Si - S N ) j (y k+1 - y k ) -j= 



Hi- 



Noting that Y/k+i) ~ ^(k) = ^(fe+i) ~ ^(fe) f° r a ^ the condition (fl3|) for stability, 
and part (i) of the Lemma, are now read from the general result about equilibrium 
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distributions of reflecting Brownian motions stated in the following lemma, and 
standard theory of Harris recurrent Markov processes ( see, e.g., [9]). 

Part (ii) is established by showing that the centered system (Xi(t) — X(t),l < 
i < N) , t > is a Harris positive- recurrent diffusion with the indicated invariant 
measure. The recurrence property has been proved in detail in [3]. The invariance 
is evident because a uniform randomization of labels relative to the centered order 
statistics is clearly invariant for the centered motion. This convergence in distribu- 
tion, combined with part (i) of Lemma [5] gives convergence in probability of Xi(t)/t 
to Sn for each 1 < i < N. Almost sure convergence can now be justified by appeal 
to an ergodic theorem for the Harris recurrent centered diffusion. □ 

The proof of Theorem [8] is completed by the following lemma, which we deduce 
from the general theory of stationary distributions for reflecting Brownian motions 
in polyhedra due to R. Williams [32] : 

Lemma 9. Let R := (Rt,t > 0) be a Brownian motion in the domain 



for some collection of K linearly independent linear Junctionals hi, with R having 
identity covariance matrix, normal reflection at the boundary, and constant drift 5 
with 



This process R has a stationary probability distribution iff a.; > for all i — 
1,...,K, in which case in the stationary state the bi(Rt) are independent expo- 
nential variables with rates 2ai, and the process in its stationary state is reversible. 

Proof. This is read from the particular case of [32J Theorem 1.2] when the matrix 
Q is identically 0. According to that theorem, R is in duality with itself relative 
to the measure p on the domain whose density function with respect to Lebesgue 
measure at x is exp(2<5 • x). Using (|16p. the linear change of variables to ?/,; = bi{x) 
shows that the p distribution of the bi(x) has joint density at (yi, . . . , yx) equal to 
c Jl^Li ex P( — 2 a iUi) f° r some c > 0. It follows that 

(17) p has finite total mass iff a% > for all i 

in which case, when p is normalized to be a probability, the p distribution of the 
bi(x) is that of independent exponential variables with rates 2a^. The "if" part of 
the conclusion is now evident. For the "only if" part we argue that if a stationary 
probability distribution p' existed, it would obviously have a strictly positive density 
on the domain. Then R sampled at time 0,1,2,... would be an irreducible Harris 
recurrent Markov process with respect p', hence p — cp' for some c > by the 
uniqueness of the invariant measure of a Harris recurrent Markov chain, and then 



{x e M. K : h(x) >0fori=l, ...K} 



K 



K 



(16) 




a, > for all i by (fl7| . 



□ 



From Theorem [8] we immediately deduce: 
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Corollary 10. For each 5 > the N-particle Atlas system with drift vector 
(5,0, ... , 0) is ergodic with average speed 5/N . The stationary distribution of 

(X (j+1) -X {j) ,1<j <N-1) 

is that of independent exponentials (Q, 1 < j < N — 1), where the rate of Q is 
26(1- j/N). 

3. The infinite Atlas model 

The infinite Atlas model can be described loosely as a countable collection of 
linear Brownian motions such that at every time point the minimum Brownian 
motion is given a positive drift of 5 > 0, and the rest are left untouched. This is an 
example of ([T]) where I = N, Si = 6 and all other 5i — 0. Throughout this section 
we take 5=1, since for our purpose here, the general case follows from the case 
when 5 = 1 by scaling. 

It can be constructed rigorously in the weak solution framework in the following 
way. Start with the canonical setup (as in Lemma [6]) of Brownian path-space: 

(18) ( C[0,oo], {^ t }o< t <oo, ), 

where {J-t} is the right continuous filtration generated by the coordinates which 
satisfy the usual conditions, and W x is the law of the Brownian motion starting 
from x £ K. We now look at the countable product of a sequence of spaces like 

OO OO 

(19) f! = C[0,oo] N , Gt=(&Tt(i), P* = (g)W% 

l l 

where the natural coordinate mapping is a countable collection of independent 
Brownian motions under P x starting from the sequence x = (x\,xi, ■ ■ ■)■ 

Let x be a sequence such that x^ > — oo and let X = (X%, X2, ■ ■ ■) denote the 
sequence of infinite independent Brownian motions starting from x. We have the 
following lemma whose proof will follow later. 

Lemma 11. Assume that the initial sequence x is arranged in increasing order 
X\ < X2 < . . . and satisfies 

(20) liminf fo""^) 2 > , 

Then, P x -almost surely, X^(t) > —00 for all t > 0, and the process 

00 „t 

(21) W*:=£ / ^.(.KwWl^W' 

i=i 7 ° 

is a {Sty-martingale whose quadratic variation (N) t = t. The stochastic exponential 
of N, given by 

(22) D t = exp (N t - t/2) . 
is hence again a non-negative {Q t } -martingale. 

We change the measure P x by using the martingale D, i.e., define Q by 
Q x \ Gt := D t ■ P x \g t , t >0. 
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By Girsanov's theorem the probability measure Q x exists and is well-defined, and 
the coordinate process under Q x is a solution of the infinite Atlas model (Q} . Hence 
Q x is the unique law of the Atlas model starting at x. 

Our aim is the following: suppose the initial points -Xi(O) < A 2 (0) < A 3 (0) < . . . 
are spread according to the Poisson process with rate two on the positive half-line, 
and we run the infinite Atlas model starting from these points. We shall prove that 
the product law of independent Exponential(2) is invariant under the dynamics of 
the vector process A of spacings given by 

(23) Ai(t) = X (l+1) (t) - X (t) (t), i = l,2,.... 

The proof is achieved through a scries of lemmas, the main argument being 
comparing the infinite Atlas model with the finite Atlas model and suitably passing 
to the limit. Throughout the proof the probability space is given by (JTUJ) . 

We start with the following lemma whose proof follows directly from Lemma [5] 
We will find it convenient to use the following notation for the operator which sorts 
a given finite vector. If a; € M. n , 1 < n < oo, define 

(24) S(x) = (x {1) ,x {2) ,...,x in) ). 

Lemma 12. Every uj € fl comprises of a sequence of processes w{t) = (t) , ui2 (t) , ■ ■ ■)• 
For any N G N, let us denote the ordered values of the processes with the first N 
indices by 

Z n (oj) = (Z? , Z?, Z»)(u) = 5((wi, wa, . . .,u N )). 

Then 

(25) D N (t) = exp (j2j Q l {XM=z» (s)}dXi(s) - t/2^ 

is a Qt -martingale. We denote by Q X N the probability measure obtained by the chang- 
ing the measure P x by the martingale D^. 

Since, under every P x , each of the Brownian motions are independent, it fol- 
lows, by applying Girsanov's Theorem, that, under Q^, the first N coordinates 
(u>i, . . . , ojn) evolve according to the finite Atlas model, while the rest of the coor- 
dinates are independent Brownian motions with the corresponding initial starting 
points. 

Let n denote the probability measure whereby the points = Ai(0) < .X^O) < 
. . . are such that the spacings X i+1 (0) — Xi(0) are iid Exponential(2). 

Lemma 13. For /j,- almost every x, the measure Q x exists and we can define 

Q-fi = J Q x d[i(x). 

Proof. The proof follows from Lemma [Til and the law of large numbers. 

We now have our main theorem in this section which proves that n is an invariant 
measure for the spacings of the infinite Atlas model. 

Theorem 14. For any K € N, and any function F : M. K — > R, which is smooth 
and has compact support, and for any time t, we have 

(26) EQ-» [F(A 1 ,A 2 , A K )(t)} = E M [F(A U A 2 , A K )(0)} . 
Here, as defined in (|23|) . is the ith spacing -XVj+i) — -^(i)- 
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Proof. [Proof of Theorem[T] Since the previous result holds for a class of functions 
which determines the marginal distributions of a sequence valued process, Theorem 
[T] follows readily for 5 = 1. The theorem for the infinite Atlas model with a general 
5 > follows by scaling. 

Proof. [Proof of Lemma fTT] Fix an arbitrary time T > 0, and consider x as in the 
assumption of the lemma. We have the following claim. 



Claim 15, 

px 



(^uj : 3 M (w) s.t. V n > M, q mf T (X n (s) - X (1) (s)) > \ 



• We prove the above claim by establishing the following 



(27) P x I lu : 3 M(lo) s.t. V n > M , q mf T (X n (s) - Xi{s)) > 
This follows by defining the events 



(28) A i :=| o m| T (X i ( S )-X 1 ( S ))<0 |, i = l,2,... 

Since X, — X\ is a Brownian motion starting from (xi — x±), which for sufficiently 
large i is strictly positive (by (|20|) ). by Bernstein's inequality (2SJ page 153], one 
can easily estimate P x {Ai) < exp(~(xi — xi) 2 /2T) for all large enough i. Again by 
assumption (f20|) . we get 

limsup(P :^ (A n )) 1/,l < limsupexp(-(a;„ - x 1 ) 2 /2Tn) 

n — >-oo n 

( 1 T • t ( X n~Xl) 2 \ ^ 

= ex P -TTTF hmmf < 1. 

\ 2J n n / 

Thus, by Cauchy's root test, it follows that the series sum P x {Ai) < oo. One 

can now apply the Borel-Cantelli lemma to obtain that P a: (limsup i Ai) — 0, which 
proves ([2T]) . and hence the required claim. 

We now prove that the process N exists and is a martingale in time interval 
[0, T]. Since T is arbitrary, this proves Lemma [TT1 Define the finite approximations 

B k(t) = Y) / l{x iW =x (1) ( s )}^i(s), 0<t<T, fc = l,2,.... 
l =i- y o 

Each Bk is a stochastic integral with bounded, progressively measurable integrands. 
It is clear that they are martingales with quadratic variations 

A- 



Jo 



-i 

= L 



{ < s < t : mm X t (s) =X (1) (s) | 



where L refers to the Lebesgue measure on the line. We shall show that the sequence 
of martingales {B^} is a cauchy sequence in the H 2 -norm , and hence has a limit 
which is denoted by N, as in (|2"Tj) . 

To see this, observe that for any n, k £ N, we get 

(B n+k ~B n ) T = h\ 0<s<T: min J*Q(s) = X {1) (s) } . 

| n+l<2<n-i-A; J 
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By Claim[l5l we see that P x { lim„ i fc^ 00 (S n +fe — -B n )r = ) = 1. It is also clear that 
(B n +k — B n )T < T, for all n, k € N. It follows from the Dominated Convergence 
Theorem that 

lim E pI ( (B n+k - B n ) T ) = 0. 

n,k — >oo 

This, by definition, shows that the sequence of martingales {B n } is a cauchy se- 
quence in the H 2 norm. Since the space of continuous martingales under that 
norm is complete, the sequence {B n } converges to a limiting martingale which we 
denote N. It also follows that (N) t = lim n _ >00 (S n )t = t. Thus, N is actually a 
{CJ t }-Brownian motion. This completes the proof of Lemma fTTl 
Remark. The condition in ([20]) is clearly loose. In fact, if we consider Ai as in 
(|28|) . all we require is that, for every T > 0, one should have 

oo oo 

J2 PX ( A ^ ^ E ex P - ^i) 2 /2T) < oo, 

i=l i=l 

which is a much weaker condition than required by (|20[) . 

The proof of Theorem [14] relies on probability estimates proved in the next three 
Lemmas [TT1 HH and [TH The second one is actually a generalization of the first. 
But, we choose to do the first separately since it is simpler and more transparent. 
But first we will do some basic computations. 

Lemma 16. IfY~ Gamma(r, A) for some r > 1, then, for any t > 0, we have 




For z > y/iX, since r > 1, one has 



(Vtz-txy- 1 < (Viz) 1 *- 1 . 
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Thus one can bound the (|29jl by 

e- Y2/2t ) < -FTre^/ 2 / (Viz)- V^ 2 ^ 
/ 1 (r Lff\ 



E 



r(r) 

r(r) iviA 

f\2-i\r/2 poo 

i^— e tA " /2 / (2 W )( r - 1 )/ 2 e- u '(2 U ;)- 1 / 2 d U ;, W - z 2 /2 

(A 2 t) r/2 etA V2 2r /2-l ^r/a-lg-™^ 



r(r) y tA 2 /2 

(2A 2 tW 2 tX 2 /2 , , . 



2T(r) 

The final identity in the lemma is due to the duplication formula: 

r(r/2) _ 2 1 -" 
r(r) V r((r + l)/2)' 

Lemma 17. For any t > 0, and a/Z positive integers N satisfying N + 1 > 16ei, 
under P ■ fi, we have the following estimate of the probability of the event that 
during time [0,t], the globally lowest ranked process is in fact the lowest ranked 
process among the processes with the first N indices. 

N 



P ■ n {X {1) {s) = Z»(s), 0<s<t}>l- CV* ( 



where C\ is a positive constant. To remind the reader, the processes [Z¥ , 1 < k < 
N) have been defined in Lemma WM 

Proof. Note that the complement of the event has the following upper bound. 

l-P- l i{X (1) (s) = Zf(s), 0<s<t} 

< P ■ n ( Ui>jv+i{ Xi(s) < Z?(s), for some s £ [0, t] } ) 
(30) <P-fi( U:>w+i{ Xi(s) < for some s £ [0,t] } ) 

OO 

< P-v(Xi(s) < Xi(s), for somes e [0,t] ). 

i=N+l 

The final bound above is the so-called union bound. 

Now we will use exponential bounds for Brownian supremums to estimate P ■ 
\x ( Xi(s) < Xi(s), for some s £ [0, t] ) which is the same as P-fi ( inf < s <t(^i(s) — Xi(s)) < 
Note that under P ■ /i, the process X% — X\ is a Brownian motion whose initial dis- 
tribution is the law of Xi(0) — X\(0). This law is Gammafi — 1,2), since it is the 
sum of iid Exp(2). Conditional on X(0) = x, such that {Xi(0)-Xi(0) = y}, 

it follows from Berstein's inequality that 

inf {Xi(s) - X 1 (s)) < ) < exp(-t/ 2 /2t). 

0<s<t J 
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Thus, for i > 2, we define Y = X;(0) — X\(0) ~ Gamma(i — 1,2), and we use 
Lemma IT6l to get 



x(i-l)/2_ e 



P ■ n\ inf (Xi{s) - Xi{s)) < < E(exp(-Y72t)) 

V Q<s<t J 

( 31 ) _ 

< v/?F(2t) 



r(i/2) 

Plugging the estimate in (|3"0|) , we get 

(32) l-P.M{^i)=^(-),0<-<t}<V5Fe« £ { -%j2T- 

By Stirling's approximation, 3 C, a positive constant, such that 

(33) T(z) > C^e-'z*- 1 ' 2 , VzeR+. 
Thus, from ([32]) we get 

^ f2^( i - 1 )/ 2 e i / 2 
l-P.M{*(i)=^0<.<f}<CV5Fe* £ rL(i-i)/ a 

t=JV+l ' ' 



i=N v 7 t=W \ / 

< 2C^e 2t+1 ' 2 (Jj^yJ > when yi^vTT) < 1/2. 
This proves the estimate. 

Lemma 18. For any any i > and all positive integers 1 < k < N € N satisfying 
N — fc + 2 > 16ei, under P ■ [i, we have the following estimate of the probability of 
the event that in the time interval [0, t], the globally lowest k ranked processes are 
in fact the lowest k ranked processes among the ones with the first N indices. 

P ■ |i {(*(!), . • ■ , X m ) (») = (Z™, . . . , ZSf) (s), < « < (} 

<34 ' >.-ai^-ir(^)" W , 

/or some positive constant Ci depending on k. For k = 1, we get back the bound in 
Lemma \ll\ 

Proof. As in the previous lemma, we bound the probability of complement of 
the event. 

1 - P ■ fi {(*(!), ■ • -,*(*)) («) = W •• • , Z%) (s), 0<s<t} 
< P ■ n ( U 4 > w+1 { < (s), for some s € [0, t] }) 

( 35 ) 

^ £ p -A*{ ^i(s) < -^f(s), for some s € [0, t] } . 

i=N+l 

Now, if [N] denotes the set {1,2,..., N}, let us note that 

= max min X\. 

h <i 2 <...<«-! ie[iV]\{ii,i2,--.,ife-i} 
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Thus, to get an upper bound on P ■ fi { Xi(s) < Z^(s), for some s£ [0, t] }, one 
can apply the union bound one more time to obtain 

(36) P ■ fj, { Xi(s) < Zg(s), for some s G [0, t] } < 

P-a[Xi(s)< min Xj(s), for some s e [0, t] |< 

»1<. ..<**_! 

some s e [0, t],i* = min{[iV]\{ii, . . . , ik-i}} ) 

il<...<ifc-l 

Now, we can count the frequency with which i* takes its possible values as the 
choice ii < t2 < ■ ■ ■ < ik-i varies in {1, 2, . . . , N}. Let g{i) be the number of ways 
to pick i\ < %2 < ■ ■ ■ < ik-i such that i* = i. Then, it is straightforward to see 

N-l 
k-1 



.9(1) = #{ {i} : h > 1 } 
, 9 (2) = #{ {z} : t! = 1, i 2 > 2 } 



AT- 2 
fc- 2 



.9(0 = #{ {*} : *i = 1,*2 = 2,.. .,»,_! = I - 1, t, > I } = /j , / A'. 

3(0 = 0, V/>fc. 
Thus, by (36]), we get 



(37) V P-/i(Xi(s)< min XAs), for some s € [0, t] 

. , „. V i6[JV]\{u,-,»fc-i} 



n<»2<---<ife 

A' 



< ^ ( ^ _ / ) p ■ A* P^( s ) < ^( s ), for some s £ [0, i] ) . 
(=i ^ 7 

Now, again as we did when deriving the last key estimate, Xi — Xi for i > I is a 
Brownian motion under P ■ [i with the initial distribution Gamma(i — 1,2). Thus, 
by Lemma [T6l we can bound 

(2t) {i ~ l)/2 e 2i 
P • /U (Xj(a) < J*Q(s), for some s G [0,f] ) < Vr7- v 



T((i-i + l)/2) 

Combining the above inequality with estimates (|33|l. (|56")) . and (13"7) . we get 
1 - P • /x { . . . , X (fe) ) (a) = (Z* . . . , Z%) (s), 0<s<t} 

^ L 2^*_ijv^ r((i-z + i)/2)- 

i=N+l 1=1 v 7 u 'III 

We will again give a loose but good enough upper bound for the infinite sum. But, 
first we need to note that, for any n > k, we have 
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Thus, it follows that 



We can use this to simplify 

1 -P- »{(X {1) , . . .,X (k) ) (s) = (Z?, . . .,Z») (s), 0<s<t} 

* 2. Z.( fc -iJ^ e r((i-z + i)/2) 

We again use Stirling's approximation to get 



i=N+l 1=1 



r((i-i + i)/2) 



i=AT+l Z=l 

oo A; 



(39) ^U-iJ^ e ^itT^-^ 1 )/ 2 ^ 072 



A- 



i=iV+i ;=i 



Since i-?>iV + l- fcfor?<fc<iV + l<i, it is clear that Aet/(i - I + 1) < 
4ei/(7V + 2 - jfe). And hence 



i — I u / . \ i—l 



1=1 \ / z=i 



By our assumption 



Aet 1 

(40) ^N=k + 2*2 



Thus 



— Z / . \ i — k 



Aet \ Aet \ , , 

< \ ^ ;— r , l<k<l, 



N-k + 2 J ~ \ \l N - k + 2 

and consequently 



z=i 



£— Z / . \ i— k 



Ei I Aet \ , Aet 
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Plugging this bound in (|39|) we get 

oo k / _ T 

e El:, 1 



k -'N-l\ ^ x (2()"-" /2 
7re 



(41) 



r((i-i + l)/2) 

i=N+l 1 = 1 v ' u " ' 

v 7 i=jV+l \ / 

N+l-k 



iV-/c + 2 / ^ \ V iV-£; + 2 

\ N-k+l 

4et \ 



by m, 



N-k+2 I 

for some positive constant Ci depending on k. This proves the lemma. 
Lemma 19. For any N € N, define fi^ to be the law under which 

X 1 (0) = 0, X i+1 (0)-X t (0)~Exp(2(l-i/N)), i = 1, 2, . . . , N - 1, 
1421 and X i+ i(0) -Xi(0) ~ £Srp(2), i = iV, N + 1, . . . , 

and all these spacings are independent. 

For any t > 0, and any three integers k < J < N satisfying 

(43) J - k + 2 > 16ei, 

we /lave i/ie following bound on the probability that under P ■ /ijv, during the time 
interval [0,t], the lowest k ranked processes among the processes with the first N 
indices are in fact the lowest k among those with the first J indices. 

(44) 

P- MAr { (Z(,...,Zi)(s) = (Z?,...,Z»)(s), 0< s <t} 



J-k+l 



4et 



N-J' 



J-k+2 J 

Proof. We follow the same line of argument as in the last two estimates. Thus 
i-p-t* N { ■ • • . Z() (s) = (zf , . . . , Z?) («), < s < t } 

N 

< P ■ Mjv ( Xi[s) < Z'l (s), for some s G [0,t] ) . 
i=J+l 

Following similar counting arguments as in (|3C5[) and (|37[) . we can bound 
l-P- MiV {(^,...,^)( S ) = (^,...,^)( S ), 0< s <t} 

(45) * /j _ A 

< ^ H(fe_J p -^A r ( < ^ some s g [0, t] ) . 
i=j+i ;=i ^ 7 

Now, under /lijv, the gap Xj(0)— JQ(Q) = ^jj where the Yj's are independent 

and Yj is distributed as Exp(2(l — j/N)). Thus, by the exponential bound used 
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before, we get 



(46) 



P ■ jt/jv ( Xi(s) < Xi(s), for some s e [0,t] ) 



< E 



exp 



V 



J= l 



Now, since each Yj is Exponential(2(l — j/N)), the random variables 

y; = {\-j/n)y 3 



are iid Exp(2), and each Y* < Yj since j < N. Thus, J2]=i Y j < Y,)~i Y j> and 
hence 



E 



exp 



V 



"5 ^ 



< E 



exp 



But, each Y* is Exp(2), and hence 



E 











exp 


;-#•)]; 


= E 


' e -Y 2 /2t 











where Y is a Gamma(z — i, 2) random variable. Thus, plugging in the bound from 
Lemma IT6l into (j46|). we derive 



P-Hn( Xi(s) < Xi(s), for some s € [0, i] ) < 



e 2 *(2*)t*- , )/2 

r((i-/ + i)/2)- 



Now, we follow approximations similar to (|39[) and (|4"Tj) for the right-hand-side on 
(|45|) to obtain 

i-p-n N {(z(,...,z£)(8) = (z?,...,zF)(8), o< s <t} 

_ e 2t (2t) (i-0/2 



JV fc 

< E E 

i=j+i i=i 



J -I 

k-1 



m-i + i)/2) 

N k / I 



2t+l/2 



< C 



J -I 

k-1 

2/ / 



'ire 



24+1/2 



< C 2 (k)e 2t (J- 1) 



E E 

i=j+i i=i 

E *(< 



4et 
J-k + 2 



i—k 



(J-k + 2) ) 

/ j=0 



by®, 



4ei 



J-fc + 2 



Now, if we call r = y/4et/(J - A; + 2), then r < 1/2 by asusmption ((43]). The finite 
geometric sum can be easily bounded as 



N-J-l 

E 

3=0 



1—7' 



N-J 



<2(l-r N - J ) 



1 - r 
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Thus, suitably altering the constant C 2 , we get our desired bound 
l-P-» fr {(Zi,...,Z£)(s) = (Z> T ,...,Z?)(s), 0<s<t} 



<C 2 {k)e 2t {J -if' 1 



This proves the lemma. 



J-k+l 



J-k + 2 



1 - 



4ei 



N-J' 



J-k + 2 



Proof. [Proof of Theorem [W] For every K < N E N, let 

Y i lf = Z$. 1 -Z?, l<i<N. 

From Corollary [TOl we know that the law of the first (N — 1) spacings under /z w 
(defined in Lemma [T9|) is exactly the stationary distribution of spacings for the finite 
Atlas model with N particles. Since, under Qn (see Lemma [12]|. the dynamics of 
the processes {X\(t), . . . , Xpf(t)} is that of a finite Atlas model independent of the 
rest of the Brownian motions, it follows by stationarity that for any t > 0, we have 

(47) e Q«-m« [F(YF, . . . , Yg)(t)] = W N [FiYf, . . . , F#)(0)] . 

We will show that for fixed t and K, as N tends to infinity, the two sides of the 
above equation converges to the corresponding sides of (I26p . This will prove the 
theorem. 

However, to do this we will need an intermediary stage where for an integer 
J < N, the dynamics of the process is according to Qj while the initial distribution 
of the spacings is either /1 or /ijy Define, for K < J < N, the following quantities: 

a = E Q N ; N [ F(X N t . }Y »)(tj\ - E«'-"* [F(Y/, . . .,Y£)(t)) , 

b = E Q,-M« . . . ; Y £){t)] - E<*'-» . . . , Y£){t)] , 

c = E^[F(Y 1 J ,...,Y^)(t)] -E^[F(A 1 ,...,A K )(t)}, and 

d = W[F(A 1 ,...,A K )(0)}-W[F(Y 1 N ,...,Yg)(0)]. 

It is clear that a, b, c, and d, all depend on t, K, J, and N, although we choose 
to suppress this dependence in the notation. Also, it follows plainly from their 
definitions combined with equality (I47|) that 

I [F(A l5 . . . , A K )(t)} - W [F(A 1 , . . . , A K )(0)] | < \a\ + \b\ + |c| + |d| . 

We will now show that a,b,c,d all go to zero if we select a sequence of J and N 
such that J, N and Nj J 2 go to infinity. This will prove the theorem. 

Step 1. [Estimate of a] For x\ < X2 < • • ■ < xk+i, define 

G(xi,x 2 , ■ ■ .,x K +i) ■= F{x 2 - xi,x 3 -x 2 ,.. - ,x K +i - x K ). 

Then, clearly, G is also a continuous bounded function. Assume that sup^ |G(x)| < 
a. Now, by the changing the measures from Qn and Qj to P, we get 

a = E Q ^~ [G(Z», Z$ +l )(t)] - E«'-"* [G{Z(, . . . , Z J K+1 ){t)] 

= E P '^ [D N (t)G(Zf, Z» +1 )(t) Dj(t)G(Zi, . . . , Z J K+1 )(t)] , 



where Dn and Dj are the Radon-Nikodym derivative processes defined in ([25 
Define the event 

r := { (Zf , . . . , Z» +l ) (s) = (Z?,..., Z J K+X ) («),V.€ [0, *] } . 
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If u) € r, by definition, P ■ /ijv almost surely Djv(t,w) = Dj(t,uj). Thus (|48j) can 
be written as 



a = E p -^[(D N (t)G(Z?,...,Z% +1 ) 



Dj{t)G{Z{,...,Z J K+l )) Ire] 

Thus, |o| is bounded above by 

E p '^ | D N (t)G(-)l r .\ + E p -^ | Dj(t)G(-)l r c | 

(49) < a (E p ^« | D N (t)l rc \ +E p -^ \ Dj(t)l r , \) 

<a(\\D N (t)\\ N + \\Dj(t)\\ N ) y/P.n N (T°)- 

The final inequality is due to the Cauchy-Schwartz inequality, where the norm \\-\\ N 
refers to the L 2 norm under the measure P ■ /i^r ■ 

Now, under P ■ fi^, both the Radon-Nikodym derivatives D]y(t) and Dj(t) are 
equal in law to exp(i? t —t/2), where B is a standard Brownian motion. Thus, it is 
straightforward to see 

\\D N {t)\\ N = \\Dj{t)\\ N =e^{t/2). 

Also, by Lemma HUl (put k = K + 1 in the statement), for large enough J and N 
we have P ■ (r c ) is less than 



C 2 (K + l)e 2t {J -1) 



K 



J-K 



J-K+l 



1 



4et 



J-K+l 



N-J' 



If we plug everything back to (j4T)|) , we see that \a\ goes to zero as J, N, and N/J 2 
go to infinity, while keeping t and K fixed. 

Step 2. [Estimate of b] Under Qj, the vector (Yj , . . . ,Y^) depends only on the 
first J processes X%, X 2 , ■ ■ ■ , Xj, and independent of Xi, i > J. Thus 

E«5 [F(Y{ J ,...,Y^)] =H(xi,x 2 ,...,xj), 

where H is bounded function since F is bounded. Thus, we have the following 
equality 

|6| = | (H (X^O), . . . , Xj(0)j) - E^( J ff(X 1 (0), . . . , Xj(0))) \ 

If and v denote the law of the vector (Xi(0), . . . , Xj(0)) under /x^r and [i 
respectively, then, since we have assume G (or equivalently F) to be bounded 
in absolute value by a, we have |6| < a \\vn — HI w Here, refers to the total 

variation norm. We will now show that — HItv goes to zero as N, J, and 
Nj J 2 go to infinity. 

Under fi^, Xi(0) = and each initial spacing Xi + i(Q) — -2Q(0), 1 < i < N — 1, 
follows independent Exp(2(l — i/N)). While, under /i, Xi(Q) — and the spacings 
follow iid Exp(2). Now the law of the vector (JQ(0), 1 < i < J) is determined 
by the first J spacings (X i+1 (0) — Xi(0), 1 < i < J) which gives us the following 
inequality 

f[ 2(1 i/N)e- 2 ^- l ' N ^ - 2 J e- 2 ^L* Xi ci.r 



\"N 



I TV 



< 



i=l 
J 



H(l-i/N)e 



2iXi/N 



- 1 



i=l 



2 J e~ 2 ^ Xi dx. 
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By an application of Cauchy-Schwartz inequality, we get 

2 

i=l 



\vn - v\\tv 



2 < 



Y[(l - i/N)e 2ixi ' N - 1 
i=l 

(50) =2 J fj(l-z/iV) 2 / cxp l-2^{l-2i/N)x\dx 

i=l "^ RJ I i=l J 

- 2 J+1 TT(1 - z/JV) / exp \-2S^{l- i/N)x\dx + \ 

7=i J ^ J { tt J 

By the standard identity J e~ Xx dx = A -1 for A > 0, we get 

J J 

»=i »=i 

.7 .7 

( 51 ) - 2 j+i jj(i - »/jv) jj^C 1 - wr* + 1 

i=l »=1 

_ (1 - l/iV) 2 (l - 2/N) 2 ... (1 - J/AQ 2 _ 
~~ (1 - 2/N)(l - A/N) ... (1 - 2J/N) ~ 
The following inequality is straightforward to prove 

e~ 2x <l-x<e~ x , for all < x < ~ log 2. 

By our assumption J, N, and N/J 2 are going to infinity, we can assume for all 
sufficiently large values of J and N that 2 J/yY < log2/2. By the previous inequality, 
we get 

e~ 2i ' N < (1 - i/N) < e-V", 1 < i < 2 J/N, 

and consequently 

(1 - 1/tV) 2 (1 - 2/AQ 2 ... (1 - J/N) 2 exp (-2g^i/jv) 
(1 - 2/JV)(l - 4/iV) ... (1 - 2J/JV) " exp (_ 2 ^ =1 2i/Jv) 



exp(~ 4 EtiV^) ^ (1 - 1/AQ 2 (1 - 2/iV) 2 ... (1 - J/N) 2 
exp(-Eti2*/iv) 



(l-2/N)(l-4/N)...(l-2J/N) ' 

Thus, we can give upper and lower bounds 

e -j(J+i)/N < (l-l/AQ 2 (l-2/AQ 2 ...(l- J/AQ 2 J(J+1)/JV 
~ (1-2/jV)(1-4/jV)...(1-2J/AT) ~ 

Combining this with (|5Tj) we see that as N, J, and N/J 2 go to infinity, we clearly 
get that |6| converges to zero. 

Step 3. [Estimate of c] This is similar to the methods we use to estimate a. As in 
Step 1, we first employ a change of measure to obtain 

c = E p -»[Dj(t)G(Zi, . . . , Z J K+l ){t)] - V p ^[D(t)G(X (1)l . . . , X {K+1) )(t)]. 

Consider the event 

(53) T(i, J, K) = { . . . , X {K+1) ) (a) = {Z(, Z J K+1 ) (s), 0<s<t} 
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On T(t, J, K), the random processes G{Z{ 1 . . . , Z , j <+1 ){s) are identical to G{X^ 1 . . . , X^ K+ i^)(s) 
in time s £ [0, t). Also, the processes Z( and Xn\ are the same in the time interval 
[0, t]. Thus the process D(t) and Dj(t) are also the same. Since, by our assumption, 
\G\ is bounded by a, the following upper bound on |c| holds 

\c\ < aE p » [(\D(t)\ + \Dj(t)\) l r .(t,j,jf)] • 

< a (\\D t \\ + \\Dj(t)\\) s/P ■ p(T°(t, J, K)), 

where, we denote by ||-|| the L 2 norm under the measure P ■ /i. 

Now, by Lemma [TBI (for k = K + 1 and N = J in the statement), for large 
enough J we get 

/ , \ J-K 

(55) P ■ fi( T c (t, J, K) )<C 2 {K+l)e 2t (J- 




Now, since t and K are fixed, as J tends to infinity, P ■ fi( T c (t, J, K) ) goes to zero. 
Finally, as in the estimate of a in Step 1, note that, by (|22]) and (|25|) . we can assert 
that 

D(t) = exp(X t - t/2), and D N (t) = exp(F t - t/2), 

where X and Y are Brownian motions. Thus ||-D(f)|| — e*^ 2 — \\DN(t)\\. If we plus 
back everything in (|54p . we get |c| goes to zero as J tends to infinity. 

Step 4. [Estimate of d.] At time zero the indices are arranged in the increasing 
order, and hence (Ai, A2, . . . , Ak) is obviously equal to (Y^ , , . . . , Y^), Since 
\F\ < a, it follows that \d\ is bounded by a times the total variation distance 
between the law of (Ai, A2, . . . , Ak ) under n and /in- A computation exactly 
similar to ([50)) gives us 

\d\\ 2 < {l~l/N) 2 (\-2/N) 2 ...{I- K/Nf 1 



{l-2/N)(l-A/N)...(l-2K/N) 

But, since K is fixed, and N grows to infinity, the right-hand-side above goes to 
zero by a logic similar to (|52|) . This proves the estimate. 

Thus, combining Steps 1, 2, 3 and 4, we have proved the theorem. 
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